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A New Method of the Corotation Radius Evaluation in our Galaxy 
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ABSTRACT 

We propose a new method for determination of the rotation velocity of the galactic spiral den- 
sity waves, correspondingly, the corotation radius, rc, in our Galaxy by means of statistical 
analysis of radial oxygen distribution in the galactic disc derived over Cepheids. The coro- 
tation resonance happens to be located at rc ~ 7.0 — 7.6 kpc, depending on the rate of gas 
infall on to the galactic disc, the statistical error being ~ 0.3 — 0.4 kpc. Simultaneously, the 
constant for the rate of oxygen synthesis in the galactic disc was determined. 

We also argue in favour of a very short time-scale formation of the galactic disc, namely: 
tf ~ 2 Gyr. This scenario enables to solve the problem of the lack of intergalactic gas infall. 
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1 INTRODUCTION 

It is recognized that the corotation resonance where the rotation ve- 
locities of galactic matter and density waves, responsible for spiral 
arms coincide, plays an important role in the galactic evolution and 
the value of the corotation radius is one of the fundamental galac- 
tic parameters. In literature one can find several methods of this 
quantity deriving in our Galaxy, for instance, by means of analysis 
of arms geometry, stellar or open clusters kinematics, HI and HII 
emissions (e.g., Lin et al. 1969; Marochnik et al. 1972; Mishurov 
& Zenina 1999; Lepine et al. 2001; Fernandez et al. 2001; Dias & 
Lepine 2005, etc.). 

In the present paper, we propose a new approach to evalua- 
tion of the corotation radius in our Galaxy. The method is based on 
the statistical analysis of oxygen radial distribution in the galactic 
disc (notice that Martin & Roy 1995 and Scarano et al. 2010 also 
mentioned the importance of the corotation effects in analysis of 
abundance gradient in external galaxies). Oxygen was used since it 
is mainly produced by SNe II which are strongly concentrated in 
spiral arms. Hence it is the most pure indicator of spiral arms in- 
fluence on the formation of radial abundance pattern in the galactic 
disc. 

As an observed material, we use the data on oxygen distribu- 
tion derived by Andrievsky et al. (2002 a,b,c) and Luck et al. (2003, 
2006) over Cepheids. Being bright and very young objects with 
precise distances, these stars give reliable information about abun- 
dances of heavy elements, close to the one in interstellar medium, 
in the significant part of the galactic disc. 

The main finding of the above papers is that the radial distri- 
bution of metallicity in the galactic disc is bimodal, i.e. there is a 
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rather steep gradient in the inner part of the disc at 5 ^ r ^ 7 kpc 
and a plateau-like distribution for r > 7 kpc and up to about 10 kpc 
(the solar galactocentric distance ro = 7.9 kpc). Hence, there is a 
bending in the slope of the distribution at r ~ 7 kpc. 

The above fine structure in the radial abundance distribution is 
very important. It indicates that in the galactic disc the distribution 
formation is caused by some non-trivial process. Mishurov et al. 
(2002), Acharova et al. (2005; 2010) developed a theory of spiral 
arms influence on the radial distribution of oxygen in the galactic 
disc. They show that the bending in the slope of oxygen distribu- 
tion is associated with the corotation resonance. However, in our 
previous papers we fitted the theory to the observations "by eyes". 
Now we propose a statistical method for the deriving of the corota- 
tion resonance location by means of analysis of oxygen distribution 
along the galactic radius. Simultaneously the so-called "constant 
for the rate of oxygen synthesis" is estimated. 



2 STATISTICAL METHOD 

Let [X/H]° b (r) be an observed distribution of any X element 
along the galactic radius r (as usual [X/H] = log(Nx /Nh)s — 
log(Nx /Nh)q, where Nx,h is the number of X element or hy- 
drogen atoms in the object, the first item refers to a star located 
at the distance r, the second one - to the Sun). On the other hand, 
let us assume that [X/H] th (r) is a theoretical distribution of the 
corresponding element which depends on some vector A, the co- 
ordinates of the vector being the sought-for free parameters of our 
theory. To fit the observations we minimize the variance a 2 over A: 
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here the summation is taken over all i — th points of galactocentric 
radius where the abundances of the element were measured, n is 
the number of observational data, p is the number of the sought-for 
free parameters. 

Unfortunately, [X/H] th (r) cannot be represented by any an- 
alytical or, at least, approximating formula. The theoretical abun- 
dance is derived as a numerical solution of a system of equations, 
the solution being dependent on a particular set of the free param- 
eters. In our approach, there are 2 target free parameters (p = 2). 
To find them the following method is used: we solve numerically 
the equations of oxygen synthesis varying the coordinates of A - 
vector within some regions and compute a net of cr(A). After that 
we construct the surface a as a function of the two free parameters 
and find the minimum (cx m = min a) which determines the best 
set of the sought-for parameters. 

To estimate the free parameters errors the confidence region 
was constructed: according to Draper & Smith (1981) the confi- 
dence contour is determined by the intersection of the surface ct(A) 
with the horizontal plane at level: 

P 



-F(p,n-p,0.95)], 



(2) 



where F is Fisher's F-statistics and the 95% level of confidence 
was adopted. 



3 EQUATIONS OF OXYGEN SYNTHESIS 

The equations which describe the evolution of oxygen synthesis in 
the galactic disc are as follows (details see in Acharova et al. 2005 
and 2010): 



/ 



(m — m w )ip(t - T m )<t>(m)dm, 



(3) 



mjj 



(4) 



= ip - J mip(t - Tm)ct>(m) dm, 
3 f 

—q^- = / (m-m w )Z(t-Tm)ip(t-T m )(j>(m)dm + 

E + fZ f -Z^--- [runo - r» 9 D—) , (5) 

where /i 9 , s are surface densities correspondingly for interstellar gas 
and stars, uo is the density of oxygen in interstellar medium, Z = 
Ho/Hg is the fraction of oxygen (hence [O/H] = log(Z / 'Zq)), 
ip is the star formation rate (SFR; ip = vu x g r \ v is a normalizing 
coefficient), <p(m) is Salpeter's initial mass function with the ex- 
ponent of - 2.35 (stellar masses in are in solar units), E is the rate 
of oxygen synthesis, / is the infall rate of intergalactic gas on to 
the galactic disc, u is the radial velocity of gas within the galactic 
disc averaged over the galactic azimuth, t is time (in Gyr), r m is 
the life-time of a star of mass m on the main sequence, tul = 0.1, 
mu — 70, m w is the mass of stellar remnants (white dwarfs, neu- 
tron stars, black holes). 

The infall rate of intergalactic gas on to the galactic disc is 
described as / = Aexp(— r/ra —t/tf) with the radial scale r<j = 



3.5 kpc (Marcon-Uchida et al. 2010) which is an intermediate value 
between rd = 2.5 and 4.5 kpc usually used in galactic nuclear 
synthesis modeling (e.g., Naab & Ostriker 2006; Fu et al. 2009; 
Schonrich & Binney2009). Constants tf and A are defined below. 

We made experiments with various abundances of the infall 
gas from Zf — 0.02 Zq to 0.1 Zq. In accordance with Lacey 
& Fall (1985) our final abundances weakly depend on the exact 
value of Zf if Zf < 0.1 Zq. Below we demonstrate the results 
for Zf — 0.02 Zq which is slightly less than the mean content of 
heavy elements in halo stars (~ 0.03 Zq, Prantzos 2008). 

The enrichment rate of interstellar medium by oxygen is rep- 
resented as follows: E = nPR, where P = 2.47 is the mass (in 
solar units) of ejected oxygen per one SNII explosion (Tsujimoto 
et al. 1995), R is the rate of SNe II events 



R(r, t) = 0.9975 J ip(r, t ~ T m )4>{m) dm, 
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the factor r\ describes the influence of spiral arms on radial oxy- 
gen pattern formation and is represented by the expression: r/ = 
f3\fl(r) — fip|Q, where fi(r) is the angular rotation velocity of 
galactic disc, Qp is the angular rotation velocity of spiral density 
waves (recall that whereas the galactic matter rotates differentially, 
i.e. fi(r) is a function of r, galactic density waves rotate as a rigid 
body, Sip — const, the corotation radius rc is determined from 
equation Q(rc) = fip), P is the constant of the rate of oxygen 
synthesis, is a cut-off factor: 0(r) = 1 if r is contained between 
the inner and outer Lindblad resonances and = otherwise. 
In our computations, we use the rotation curve based on Clemens 
(1985) data adjusted for the scale r = 7.9 kpc. 

The last term in parentheses in Equation (5) describes the 
turbulent radial diffusion of heavy elements. To derive the diffu- 
sion coefficient D we use the simple gas kinetic approximation 
(Mishurov et al. 2002). 

The above system of equations may be divided into 2 groups: 
i) equations (3,4) which describe galactic disc formation, and ii) 
the "chemical" equation (5) of oxygen synthesis. The target 2 free 
parameters, tip and j3, enter only the equation (5). Hence, the first 
group of equations can be solved independently of the last equa- 
tion. Therefore, the algorithm of solving the above equations is as 
follows. 

At first step, independently of Equation (5) we solve equa- 
tions (3,4) by means of Lelevier method (Potter 1973) and find the 
time evolution of the radial profiles of gaseous and stellar densi- 
ties, u g (r, t) and u s (r, t) (simultaneously SFR, ip(r, t) is derived). 
To perform this part of modelling, we have to specify a value of 
the time-scale of gas infall on to the galactic disc, tf, and adopt 
an expression for u(r). After that, the constants v and A, entering 
these group of equations, are obtained by fitting the normalising 
conditions that at present epoch (t = Tp, — 10 Gyr - the age of 
galactic disc) and solar position stellar and gaseous densities have 
to be equal to: /x s (?"o,Td) = 40 M pc~ 2 and u g (ro,Tp>) = 10 
M pc~ 2 (Haywood et al. 1997). Taking into account that the part 
of stellar remnants is of the order of 10% of fi s (our experiments 
confirm this supposition) we have that at present time and solar 
galactocentric distance summary density has to be about 54 Mq 
pc~ 2 - a value close to Marcon-Uchida et al. 2010. However, be- 
low we also discuss other normalising conditions. 

At the second step, we solve the chemical equation (5) for a set 
of ftp and /3 (as boundary conditions we demand that the solution 
has to be finite at the galactic centre and the galactic disc edge 
which we dispose at r = Rg = 35 kpc). Repeated this procedure 
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we find the minimum of a — a m , hence estimate the best values 
for Qp and /3, then compute their errors as it was described in Sec. 
2. At last we derive the location of the corotation resonance rc- 



4 STEP 1: GALACTIC DISC FORMATION 

In what follows we examine 2 limiting cases for t / : t / = 2 Gyrs - 
rapid disc formation, and tf = 7 Gyrs - slow disc formation. For 
the radial gas flow 3 models were studied. 

1) No radial flow of gas: u = (models M20 and M70 corre- 
spondingly for the above two rates of disc formation). 

2) Gas radial inflow: case a) u = — 0.3r km s" 1 , if r < 3.3 
kpc and u = — 1 km s -1 , if r ^ 3.3 kpc (models M2a and M7a); 
case b)u = — O.lr km s _1 (models M2b and M7b); 

3) Gas radial outflow: u = +0.5 km s _1 (models M20 and 
M70). 

In Fig.l evolution of gaseous radial profiles with time is 
shown. In significant part of galactic radii current gaseous density 
distributions happen to be close to the observed ones and profiles 
derived by Portinari & Chiosi (2000), Fu et al. (2009), Naab & 
Ostriker (2006) and Schonrich & Binney (2009) . Comparing the 
distributions we cannot prefer a certain model only on the basis of 
density profile (models M20 and M70 have holes in the galactic 
centre but they do not describe oxygen radial distribution well). 

So, we need some additional arguments which enable us to 
make a choice among various scenarios of galactic disc formation. 
Such information comes from observations: according to Sancisi 
& Fraternali (2008) and Bregman (2009) current rate of total (inte- 
grated over the galactic disc) gaseous mass falling on to the galac- 
tic disc is of the order of 0.1 - 0.2 Mq yr _1 whereas the total star 
formation rate is ~ 1 - 2 Mq yr _1 (Robitaille & Whitney 2010) 
although one can meet a higher level of SFR ~ 5 Mq yr _1 (Hart- 
mann 2007). The above values may be considered as global nor- 
malisations which are complementary to the particular ones at solar 
position. 

In Table 1 we give the computed total rate of forming stars, 
M s = 2-7T J^ G iprdr, and infall gas mass, Mf = 27r J^ g frdr at 
the present epoch. From the Table one can see that only models of 
rapid galactic disc formation (tf ~ 2 Gyr) demonstrate low current 
infall rate, Mf ~ 0.1 - 0.2 M yr _1 , and SFR at level M s ~ 1.5 
Mq yr _1 , the both quantities fit the above observed values. Sce- 
narios of slow disc formation (tf ~ 7 Gyr) lead to greater current 
SFR, Ms ~ 3 - 4 M yr" 1 , but they require the infall rate of the 
order of magnitude higher than that observed by Sancisi & Frater- 
nali (2008) and Bregman (2009). Obviously, one should be prudent 
in such estimates since future observations could discover such a 
huge amount of baryonic matter falling on to the galactic disc. But 
at the moment we do not have such data, at least for our Galaxy. 

Let us look at the problem from another angle. There is a be- 
lief that star formation is mainly fuelled by the gas infall on to the 
galactic disc. Hence, the total rates of gas infall and star formation 
should be close to each other in order to the interstellar gas not to be 
consumed due to star formation. But as we wrote above, at present 

1 Strictly speaking in a full formulation of the problem we have to solve 
simultaneously gas dynamical equations which take into account the effects 
of the central bar, spiral arms, etc. However, such complicated task is still 
beyond the scope of approach usually used in theories of nuclear synthesis 
in galactic discs. Instead of that, following Lacey & Fall (1985), Portinari 
& Chiosi (2000) we examine several model representations for the radial 
velocity u. 
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Figure 1. Evolution with time of the profiles of gaseous density: dotted line 
- for t = 2; dashed - dotted line - for 4 = 4; dashed line - for t = 8; solid 
line for t = 10 Gyr. 

Table 1. Constants A (Mq pc" 2 Gyr" 1 ) and v (M~ 5 pc Gyr" 1 ); 
integrated rates of gas infall, Mf, and star formation rate, M s (Mq yr -1 ) 
and minimal u m derived in various models 
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0.14 


M f 
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0.050 
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0.060 


0.066 


0.050 


0.043 


0.086 


0.068 



epoch in our Galaxy the total star formation rate is at least 5 times 
higher than the one of gas infall. 

To explain this problem, in Fig. 2 we show evolution of the 
ratio of gas infall rate on to the galactic disc to star formation rate, 
f/tp, with time as a function of r for 4 models. From the Figure it 
is seen that, at each moment of time there is a region in the galactic 
disc where f/ih > 1. This means that the interstellar gas is not 
fully consumed in the Galaxy (from Fig. 1 one can also see that the 
gaseous density is not equal to zero at any galactocentric radius). 

Besides, as it was expected, in scenarios with gas radial inflow 
the region, where f/tp > 1, is shifted closer to the galactic centre. 
Indeed, gaseous density decreases with galactocentric radius, star 
formation rate also decreases with r. Hence, due to interstellar gas 
radial inflow the region with depressed star formation rate moves 
closer to the galactic centre and this results in the ratio f/ip in- 
crease. In addition gas radial inflow replenishes the inner part of 
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Figure 2. Evolution with time of the ratio of infall rate to SFR (line types 
correspond to the time moments as in Fig. 1). Horizontal dashed lines sep- 
arate regions in the Galaxy where / > ip. Here a reserve of gas is formed 
for the future star formation. 

the galactic disc with building material for star formation fuelling. 
So, as we can see, at each epoch of galactic evolution the reserve of 
gas forms into the galactic disc and it fuels star formation process. 

How should we change the normalisations of the current 
gaseous and stellar densities at solar distance in order to obtain 
the global star formation rate at level, say, M s ~ 5 M Q yr" 1 ? 
Our experiments give the following results. In scenario of rapid 
galactic disc formation (tf ~ 2 Gyr) and no gas radial inflow, 
gaseous and stellar densities happen to be unrealistically large: 
fi s (r ,T D ) ~ 155 Mq pc" 2 ; /i 9 (r ,T D ) ~ 30 M pc" 2 albeit 
the rate of gas infall Mf ~ 0.5 Mq yr" 1 (in case of gas radial 
inflow like in models "b" /x s (rn, To) ~60 Mq pc" 2 ; Mf ~ 0.3 
Mq yr" 1 however n g (ro, To) keeps at very high level ~ 30 Mq 
pc -2 ). In scenario of slow galactic disc formation (tf ~ 7 Gyr) 
and no gas radial inflow we have: fX s (To,Tn) = 75 Mq pc" 2 ; 
Ms( r o, T D ) = 14 M pc -2 but the rate of gas infall Mf ~ 3.1 
Mq yr" 1 is of the order of magnitude higher than the observed 
value (if we take into account gas radial inflow like in models "b" 
we derive: fj, g (r ,T D ) ~ 25 M Q pc" 2 ; fi a (r ,T D ) ~ 40 M 
pc" 2 but Mf ~ 2.6 Mq yr" 1 again is too high). Perhaps such sit- 
uation can be implemented in some galaxies but it is not supported 
by observations for our Galaxy. 

There is one additional argument in favour of rapid disc for- 
mation scenario: according to Allende Prieto (2010) and Zasov & 
Sil'chenko (2010) galactic discs were formed rapidly, at least last 8 
Gyr their sizes did not change significantly. 



5 STEP 2: ESTIMATE Q P ,R C AND /3 

Despite the conclusion made in the previous Section below we give 
the results for all the above models. So, at the next step we solve 
the chemical Equation (5), compute the variance <r 2 and find its 
minimum which determines the best values of Sip and /3 (to catch 
the minimum of the variance we varied Sip and /? in a wide region: 
10 < S1 P < GOkms" 1 kpc" 1 ; 0^/3^ 0.6 Gyr.). 

As an observed material we use the data of oxygen distribu- 
tion along the galactic radius, derived over Cepheids, which were 
averaged within bins of 0.5 kpc width (see Acharova et al. 2005; 
2010). The corresponding distribution is drawn in Fig. 3 at step 




6 8 10 6 8 10 

r (kpc) r (kpc) 



Figure 3. Comparison of theoretical and observed oxygen radial distribu- 
tions for models M20 and M7b. Solid lines are for the best values of Qp and 
P, dashed lines - for (fi P + A fip ; + A ) and (fip - An p ; j3 - Ag). 








Figure 4. An example of the surface — <r (Sip , /3) in case of M20. For better 
visualisation we draw the figure "bottom up". The best estimates for the 
parameters are: Sip = 33.2 km s" 1 kpc -1 ; f3 = 0.016 Gyr. 

0.25 kpc in the region from 5.25 up to 10 kpc (in total, the num- 
ber of observed points is n — 20). The bars in the Figure are the 
rms scatters of the mean abundance within the bin. To illustrate our 
method, in Fig. 4, we show an example of the surface a(Qp, 0) for 
model M20. To estimate the errors of the sought-for free parame- 
ters, we compute the section of the above surface by the horizontal 
plane at a = a c — 1.181o- m (see equation (2), in our case the 
value of Fisher's statistics happens to be F(2, 18, 0.95) ~ 3.55, 
Draper & Smith 1981). The section is demonstrated in Fig. 5 by 
the shaded figure for model M20. Its ellipse-like border delineates 
the confidence region. The semi-axes of the ellipse give the errors. 
The results of variance computations are shown in Table 1. Con- 
sider them in more detail. 

First of all, by means of Fisher's statistics analysis the sig- 
nificance of gas radial inflow inclusion into consideration can be 
estimated. To do that we have to compare models "a","b" and "O" 
with the ones without gas radial inflow separately for the above two 
scenarios of galactic disc formation. The result is as follows. 

i) In scenario of slow disc formation the best model is M7b. 
Other models have significantly higher values of o m and they are 
to be rejected. 

ii) In scenario of rapid disc formation two models, M20 and 
M2b, have the same values of o m = 0.050, other ones should be 
rejected. However, we cannot decide between M20 and M2b: the 
null hypothesis (u = 0) cannot be rejected. In other words, in this 
scenario inclusion of the gas inflow is insignificant. 

So, below we discuss two models - M20 and M7b. Let us 
compare them and test the null hypothesis that the correspond- 
ing variances do not differ significantly against the alternative one 
that they do. For this, according to Martin (1971) we should check 
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Figure 5. The shaded region is an intersection of the surface a(Qp , /?) from 
Fig. 4 by the horizontal plane a = cr c = 0.059. This region determines the 
confidence one at probability level 95 % (see Sec. 2). Its border gives the 
estimate of the parameters errors: = ±0.001; Aq p = ±1.4 km s 
kpc -1 . 

if the value F = irnin OMwImin 0\M7&) 2 lies in the range: 
l/F a /2 < F < F a /2 (a is the significance level and we adopt that 
the degrees of freedom are approximately the same in the both alter- 
natives). From Table 1 we derive F « 1.35, whereas for a = 0.02 
F a/2 = F(18, 18,0.99) w 3.14 (Draper & Smith 1981). Hence, 
the above inequality is satisfied and we conclude: the difference 
between models M20 and M7b is undetected statistically (on the 
basis of available observational data of oxygen radial distribution 
and a rather simple theoretical approach). At the moment, the only 
way to discriminate between them is as it was discussed in Sec. 4. 
Finally, our estimates for the target parameters are as follows: 

\)M20: rip = 33.2±1.4kms -1 kpc" 1 (the corotation radius 
r c = 7.0 ± 0.3 km s -1 kpc -1 ); fi = 0.016 ± 0.001 Gyr. 

2) M7b: ftp — 30.5 ± 1 km s -1 kpc -1 (the corotation radius 
rc = 7.6± j km s -1 kpc" 1 ); fi = 0.036 ± 0.002 Gyr. 

In Fig. 3 the best theoretical oxygen distributions for the above 
two models and the distributions for the upper and lower values of 
Q.P and fi are shown. From the figure it is seen that our theory 
explains the bimodal radial distribution of oxygen in the disc of our 
Galaxy rather well. 



6 DISCUSSION 

In the present paper, a new method for evaluation of the corota- 
tion radius in the Galaxy is developed. Our approach is based on a 
statistical analysis of the bimodal structure of oxygen radial distri- 
bution in the galactic disc determined over Cepheids. By means of 
treatment of observational data we derive that the corotation reso- 
nance happens to be situated at rc ~ 7.0 — 7.6 kpc depending on 
the rate of intergalactic gas infall on to the galactic disc (the statis- 
tical error is ~0.3 kpc). The above value for the corotation radius 
is close to the solar galactocentric distance ro = 7.9 kpc. Simulta- 
neously the constant for the rate of oxygen synthesis in the galactic 
disc fi ~ 0.016 — 0.036 Gyrs was obtained. 

We also argue in favour of a short time-scale formation of the 
galactic disc, namely: tf ~ 2 Gyrs. This scenario enables to solve 
the problem of the lack of intergalactic gas infall, i.e. the very low 
present integrated rate of gas infall on to the disc, Mf ~ 0.1 — 0.2 
Mq year -1 observed by Sancisi & Fraternali (2008) and Bregman 
(2009), whereas the integrated star formation rate is expected to be 
M s ~ 1 - 2 M year" 1 (Robitaille & Whitney 2010). Higher 
level of the current global star formation rate, M s ~ 5 Mq year - , 



needs too high rate of gas infall on to the galactic disc or unrealis- 
tically high gaseous and stellar current densities at the solar galac- 
tocentric distance. 
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